Holon Pair Condensation and Phase Diagram of High T c Cuprates 
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A possibility of holon (boson) pair condensation is explored for hole doped high T c cuprates, 
by using the U(l) slave-boson representation of the t-J Hamiltonian with the inclusion of hole-hole 
repulsion. A phase diagram of the hole doped high T c cuprates is deduced by allowing both the 
holon pairing and spinon pairing. It is shown that the spin gap size remains nearly unchanged below 
the holon pair condensation temperature. We find that the s-wave holon pairing under the condition 
of d-wave singlet pairing is preferred, thus allowing d-wave hole pairing. 



I. INTRODUCTION 



It is believed that the high T c superconductivity arises 
as a result of Bose condensation of doped holes. It is still 
an unresolved problem to predict a satisfactory phase di- 
agram of involving the bose condensation of the doped 
holes and the normal state pseudogap. The normal state 
pseudogap has been observed in various experiments 
|y [Tot : NMRjEfl, neutron scattering ||, c-axis optical 
conductivity |t[ 5|, heat capacity Q , in-plane resistivity 
and photoemission Jj|[l0| . A universal dependence of 
superconducting critical temperature T c on hole doping 
rate is observed for various high T c cuprates, by manifest- 
ing a smooth increase of T c up to the optimal doping rate 
and a decrease of T c beyond it. This observation is well 
fit by a universal relation, T c /T c max = 1 - 82.6(ir - 0.16) 2 
where T 1 ™* is the maximum critical temperature at the 
optimal doping rate of x — 0.16 p| . 

Based on the t-J model Hamiltonian, the normal state 
pseudogap is regarded as the spin gap |I^-[l8|]. Earlier 
superconductivity in the underdoped cuprates was un- 
derstood as a simultaneous presence of the spin gap and 
single boson (holon) condensation ]l3|-|l^]. Recently an 
SU(2) slave-boson theory of the t-J model Hamiltonian 
was proposed to allow symmetry at both half filling and 
finite doping In this theory jl9| single boson 

(holon) condensation can be either completely destroyed 
or reduced due to a low lying fluctuation (soft mode) 
in association with the SU(2) rotation, and thus boson 
(holon) pair condensation is suggested pcfl . Earlier var- 
ious preformed-pair scenarios |2l]-|2"3| were proposed. In 
one of the scenarios , the preformed pairs become lo- 
cally available below the pseudogap temperature, and the 
critical temperature of superconductivity is determined 
by phase ordering of the preformed pairs. On the other 
hand, numerical calculations |p^-|27[| have been made to 
study the paring of doped holes in the t-J clusters and t- 
J ladder systems. A recent study for doped t-J three-leg 
ladders |28| revealed hole pairing due to the coupling of 
a Luttinger liquid to the insulating or doped spin liquid. 
In view of great interest in the role of doped holes for 
high T c superconductivity, we investigate a possibility of 
boson pair condensation by introducing a modified t-J 



Hamiltonian and derive a phase diagram in the plane of 
temperature vs. hole doping rate. 



II. MEAN FIELD HAMILTONIAN AND FREE 
ENERGY FROM A MODIFIED t-J 
HAMILTONIAN IN THE SLAVE-BOSON 
REPRESENTATION 

A necessity of introducing the electrostatic hole-hole 
repulsion is stressed in recent numerical studies of a spin 
ladder system ^7j,^9|. In order to account for a reason- 
able hole pair binding a large value of Coulomb repulsion 
energy V between two nearest neighbor (NN) holes was 
needed, i.e., V = e 2 /(e OD c ) = O.leV |tJ with c ~ 3.8A, 
the lattice constant and f M = 30 ~ 40, the dielectric 
constant |30[| . Ignoring this repulsion energy, the bound 
state of the hole pairs was found to be excessively robust 
& 

In the local U(l) symmetry conserving slave-boson rep- 
resentation, we write the full t-J Hamiltonian of the two 
dimensional systems of antiferromagnetically correlated 
electrons by including a hole-hole repulsion term, 



H = 1 E (ftMh 



h.c 



+v^24b i b]b j -(x j2flfi° 



(i) 



with Sj = l/2fj <7 a a/ sfip. Here the local constraint of sin- 
gle occupancy, £ CT fLfi<r + b i b i = 1 is assumed. fl(f t a) 
is the spinon creation (annihilation) operator and bi(b\), 
the holon annihilation (creation) operator. The third 
term represents Coulomb repulsion between holes in the 
NN sites; according to an experiment |n| the hole- 
hole repulsion decays rapidly with distance. The second 
bracketed term above implies that the NN configuration 
of two holes is energetically more favorable than other 
possible configurations. This is because the two holes in 
the NN sites break 7 bonds compared to 8 bonds for other 
configurations. This is evident from the separate inspec- 
tion of the two attractive interaction terms, JSi ■ Sj and 
— {J/A)nirij. For the latter we write, 
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(id) 
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= ™£/U + ££6k-i£*i*M- (2) 



where the local constraint of single occupancy is con- 
sidered. The effective attraction between the NN holes 
arises from the last term of the equation above. 

Applying the Hubbard-Stratonovich transformation 



Jl7t in Eq. ([!]), both the Heisenberg term and the hop- 
ping term are converted into linearized terms involving 
the hopping order field \ji — (8t/3 Jbjbi + fj^fia) for the 
exchange interaction channel and the spinon pairing or- 
der field A£ = {fj\fi[ — fj{fi\) for the pairing channel. 
Such reduction to the two channels is made by ignor- 
ing the direct (Hartree) channel based on the assump- 
tion of paramagnetic states for each site, i.e., (S^) = 
]l7t . Thus long-range antiferromagnetic fluctuations are 
ignored. The resulting effective Hamiltonian is then, 
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where a Lagrange multiplier field Xi is introduced for the 
local constraint of single occupancy for both the spinon 
and the holon. The quartic holon term (the second term 
in Eq. (|) above), ^ £ (jJ> (&}&,)(&&.) is repulsive @. 
It is important to realize that this term corresponds to 



J 



the holon exchange interaction, but not to the direct (for- 
ward) interaction nor to the holon pairing interaction. 

Thus allowing only the holon exchange channel, we lin- 
earize the quartic holon term (the second term in Eq. (^)) 
as 
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We decompose the effective holon attractive energy term (the third term in Eq. (|J)), that is, — J2(i ~ V)b\b^bibj 
with < V < J/4 into the direct, exchange, and pairing channels [^3). By introducing the Hubbard-Stratonovich 
transformation for the resulting holon pairing term and the above linearized holon exchange term, we obtain 
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r 



with (i 



^0 — J/4 and = — 3J/4, where A^ is now 



the Hubbard-Stratonovich field for holon pairing. Here 
and /ig are the effective chemical potentials for the 
spinon and the holon respectively. Obviously the two 
scalar fields, Xji = (8t/3Jb^b t + f] a f la ) and (bjfy) in 
Eq. (||) above are not independent. We allow a linear 
relation between the spinon hopping and the holon hop- 
ping order parameters, i.e., (fj^fjer) = viH^j)- Thus we 
rewrite Xji = (8i/(3 J) + 2i])(b%) or 



8*/ (3 J) + 2?] 



Xji- 



(6) 



A self-consistent determination of the hopping ratio r\ 
will be discussed later. 



Substitution of Eq. (||) into Eq. (|J) results in the 

cancelation of terms (proportional to ^j) involving the 
holon exchange channel in Eq. (^J) . By defining the pair- 
ing order parameters A( = {fj\fj+ r ,i - fjlfj+r,l) and 
A£ = {bjb'- +r ) with r = x or y, and allowing the uniform 
hopping order parameter, Xji = X-, the resulting mean 
field Hamiltonian is 
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where 



and 



J = 



3Jr/ 2 /2 + 4r?£ + J/4-V 
(27? + 8t/(3J)) 2 



t = 
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2r; + 8i/(3J) ' 



(8) 



(9) 



From the hopping ratio, rj 



where x is 



(^) - 

the hole doping rate, we note that in the half-filling limit, 
i.e., x — > 0, J in Eq. (||) approaches 3J/8. As a result 
this satisfies the SU(2) symmetry [Q of the spinon pair- 
ing order field and the spinon hopping order field, that 
is, |A^| = x at x = 0. Following a similar consideration 
to that of Lee and Nagaosa |TJ| , we obtain the effective 
spinon mass, to/ -1 = c 2 -^x ( c o is the lattice constant) 
and the effective holon mass rrib, 



mb 
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(10) 



from the use of Eq. (@) with Eq. (|). 

By allowing d-wave pairing for spinons (A^ = — = 
A-^) and s-wave (A b = A b = A b ) or d-wave pairing 
(A b = — A b = A b ) for holons, we obtain the Fourier 
transformed mean field Hamiltonian from Eq. (Pi), 



H 



MF 



2N 



J| x | 2 + ^|A/| 2 + (^-y)A 



b\2 



NX 
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Here e 



/ 



3 J 



XI k, 



= -Ztxik, V f = Mo + and 
[i b = /Xq + A with 7^ = cosk x + cosk y . A* k = ^-A^ifk 
for d-wave spinon pairing with ip k = cosk x — cosfcy, 
A^ = 2(j — V)A b jk for s-wave holon pairing and 
A^ = 2(j — V)A b tpk for d-wave holon pairing. N is 
the total number of lattice sites. From the Bogoliubov 
transformation for both the fermion (spinon) and boson 
(holon) operators, we obtain the quasi-particle excita- 
tion energies, E k — \J (e{ — (if) 2 + A{ for spinons and 



,b\2 



for holons. The spinon pair- 
ing gap (or spin gap) is E[(k) = \A f k \ = \^-A f tp k \. We 
note that the upper limit of |A|| is |e| — fi \. For both 
the s-wave holon pairing and d-wave holon pairing, the 
maximum (£^ ax ) of the holon quasi-particle excitation 
energy E\ occurs at k = (jr, ir), whereas the minimum 
(-E^in) 01 occurs at k = (0, 0) (for the s-wave holon 
pairing ^ ax/min - y/(At X T M & ) 2 - [(J = W) A"] 2 and 
^max/min = l 4 *X T A (where [i b < 0) for the d-wave 
holon pairing (— for max. and + for min.). 

The mean field free energy at doping rate x is obtained 
to be, 



Fmf(x,A',A")/7V = 2J| x | 5 



3J . f i ,J 
— A ; + (- 
4 l 2 



2V)A b 



2T In [cosh(/3El /2) In [sinh((3E b k /2)] 



X[i' 



(12) 



The minimization of the mean field free energy with re- 
spect to the scalar fields, x, A /, and A h will be numer- 
ically made as a function of temperature T and dop- 
ing rate x |35[ |. From the computed mean field values 
of the order parameters (x, A^ and A b ) and the chem- 
ical potentials and fi b ), we determine the hopping 
ratio 77 self-consistently with the use of J2a(fLfj<?) = 
-mT, k [lk(4-^)/ Ef k]^HPE f k /2) and (&t 6j ) = 

47v E fe bk(4 H b )/E b ] coth(/3^/2). 



III. COMPUTED RESULTS OF BOSON (HOLON) 
PAIR CONDENSATION AND SPIN GAP 
TEMPERATURES; PHASE DIAGRAM OF T-x 
PLANE 



In Fig. 1 we show as a function of doping rate the 
predicted spin gap (Ef(ir, 0) = ^-A^), the rescaled hop- 
ping order parameter |l3| (x multiplied by ^) and the 
holon pairing order parameter A b respectively at various 
selected temperatures. The self-consistently determined 
result of r\ is found to be 77 ~ 0.2(1 — x)/x with x > 0.01 
for J = 0.4/j. For the present study, J = 0.4t Q and 
V = 0.24999 J were chosen to best fit the experimental 
phase diagram pTj,nilB6[ in the plane of T vs. x. The 
predicted spin gap E^(tt, 0) (y-axis on the left hand side) 
at k = (ir, 0) decreases with the hole doping rate x at all 
temperatures. We find that the minimum of the free en- 
ergy occurs only with the s-wave holon pairing, but not 
with the d-wave holon pairing. The d-wave holon pairing 
is found to be unstable. In Fig. 1 the s-wave holon pair- 
ing order parameter A b (y-axis on the right hand side) is 
displayed at each selected finite temperature. It is seen 
to increase with x at all temperatures. We now examine 
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0.05 0.1 0.15 0.2 0.25 0.3 

x 

FIG. 1. The computed spin gap Eg (on the left axis) at 
k = (tt, 0), the hopping order parameter \ multiplied by 3 J/4 
(on the left axis) and holon pairing order parameter A b (on 
the right axis) at various selected temperatures. 

the variation of the spin gap with temperature below and 
above the holon pair condensation temperature, T MF in 
the underdoped region. As an example, let us choose the 
doping rate of x — 0.1, at which T^ F is obtained to be 
0.028J. Two of the three nearly indistinguishable lines 
of the spin gap E[(tt,0) at T = 10" 10 J and T = 0.02J 
(both of which represent temperatures below Th F ) indi- 
cate nearly identical spin gap sizes. Thus the spin gap 
size remains nearly unchanged below the holon pair con- 
densation temperature. The remaining four lines which 
correspond to temperatures above T^ F show a consistent 
decrease in the spin gap size as temperature increases. 

In Fig. 2, we present the computed phase diagram in 
the T - x plane. The holon (boson) pair condensation 
temperature X" MF and the spin gap temperature T^ F are 
plotted as a function of doping rate x, including the spin 
gap Ef(w,0) at T = T MF . The mean field holon pair 
condensation temperature is found to be less than 0.1J 
in the underdoped region. In this region both the d- 
wave spinon pairing with (/^ /t.) ^ and the s-wave 
holon pairing with (bibj) =/= coexist below the mean field 
(critical) temperature (T c — T^ F ). For the hole (but not 
the holon) pairs, we have (c^Cji) = (b\bj) (fcfji) / 
in the mean field approximation. For clarity we would 
like to stress that the term, 'holon' refers to the boson 
of spin and the term, 'hole' is the fermion of spin up 
and spin down. Thus this allows for the condensation 
of the d-wave hole pairs by satisfying the s-wave holon 
pairing for (b\b^) and the o?-wave pairing for (/-jf/jj.) in 

( c »T c ji) = ( b i b ])(fiffjl)- Thus t mf at each doping rate 
corresponds to the mean field critical temperature of the 
d-wave superconducting phase transition. 

The mean field spin gap (pseudogap) temperature T^ F 
may be regarded as the 'center' of crossover region which 
arises as a result of gauge fluctuations |]l6|,^8). The pre- 
dicted pseudogap (spin gap) temperature T^p is seen to 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
x 

FIG. 2. The phase diagram of holon pairing and spinon 
pairing. The AM (anomalous metallic), SG (spin-gap), FL 
(Fermi liquid), and dSC (d-wave superconducting) phase re- 
gions are exhibited with the holon-pair condensation temper- 
ature Tmf an d the d-wave spinon pairing temperature T MF . 
The spin gap Eg with k = (tt, 0) at T = T MF and the mean 
field single holon condensation temperature denoted as Tg E 
are also displayed. The normal state gap E£ disappears at 
T MP = T MF with x cr ~ 0.19. The shaded region represents 
the region at which the holon pair condensation is unstable. 



smoothly decrease with the hole doping rate as shown 
in Fig. 2. This is consistent with the recent experiments 
|p|,^|jic[| . The spin gap (tt, 0) at the s-wave holon pair 
(d-wave hole pair) condensation temperature T^p is pre- 
dicted to show a rapid decrease with x as is shown in 
Fig. 2. This spin gap (pseudogap) vanishes at the criti- 
cal doping rate of x CI ~ 0.19, at which T^ F approaches 
Tjjj F . We find that the position of the critical doping rate 
is sensitive to the choice of V. Unlike T^ F , Tjjp is nearly 
independent of V. The critical doping rate x CT decreases 
as V gets smaller. With the choice of a vanishingly small 
value of V, we find that the value of holon pair order 
parameter A fc (and thus T MF ) is excessively large even 
at a small doping rate. This is consistent with a recent 
numerical study |||. The choice of V = 0.24999 J fits 
best the experimental value of the observed critical dop- 
ing rate of x cr — 0.19 |l^,^|. For this case the effective 
NN attractive energy is about V c g = 4- — V ~ 10~ 5 J. 
Earlier the mean field single boson condensation temper- 
ature was reported to be Tg E ~ 2irxco~ 2 / nib ]37| , |l6|] . It 
is displayed in Fig. 2 for comparison with our computed 
boson pair condensation temperature (T c = T^p). Ob- 
viously we find a wide difference between the two bose 
condensation temperatures of the holon-pair boson and 
the single holon boson. The observed optimal doping 
rate of x^ p — 0.16 [0 is known to be smaller than the 
critical doping rate of 0.19 [|l0|,^6). In disagreement with 
observations, our predicted boson pair (not single boson) 
condensation temperature above the critical doping rate 
of x cr = 0.19 becomes larger than the pseudogap (spin 
gap) temperature. Earlier it was suggested that the 
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FIG. 3. The effective holon chemical potential fJ. b (T = 0) 
(on the right axis) and the inverse of the holon compressibility 
l/«6 (on the left axis) as a function of doping rate x. 

overdoped region may correspond to the region where 
the mean-field solution is not stable with respect to the 
gauge fluctuations |3^] and that the system can probably 
be described by the Fermi liquid theory 16 35^. How- 
ever, it remains to be seen for a rigorous verification in 
the future. 

In Fig. 3 we present the effective holon chemical po- 
tential // at T = (y axis on the right hand side) and 
the inverse of the holon compressibility 1/ Kb (y axis on 
the left hand side) as a function of doping rate x. The 
computed spinon chemical potential decreases monoton- 
ically with doping rate x. On the other hand, the holon 
chemical potential decreases up to x ~ 0.07 and starts to 
increase beyond this value. Allowing that the inverse 
of two dimensional holon compressibility k& holds for 

v** = ^% = St£t 63 ( where N is the numbcr 

of bosons and A is the area), we find that 1/k{, is nega- 
tive below x ~ 0.07 and becomes positive for x > 0.07. 
This suggests that in the region of x < 0.07 the holon 
pair condensation is unstable ^3] and it is stable only 
for x > 0.07. The shaded region in Fig. 2 represents 
the region at which the holon pair condensation is unsta- 
ble. It is well known that the superconductivity arises 
for x > 0.05 JyJ. Interestingly recent numerical studies 
p6| , p8[ of t-J three- leg ladders revealed that there exist 
a critical doping rate of x ~ 0.06 beyond which the hole 
pairing increases rapidly as hole density increases. It is 
of great interest to see in the future whether there exists 
any possible relevance between the present result with 
the prediction of the t-J three-leg ladders study. 



IV. CONCLUSION 

A possibility of boson (holon) pair condensation was 
described with the inclusion of a repulsive interaction 
between the NN holes in the t-J Hamiltonian. A phase 
diagram of the hole-doped high T c cuprates is derived 
by considering both the holon pair condensation and the 



spinon pairing gap (spin gap) . In the mean field approxi- 
mation, we note that (c^cji) = (&J^)(/it/j'i) f° r * ne n0 ^ e 
pair. To avoid confusion, we would like to point out that 
the term, 'hole' stands for the fermion of positive charge 
+e with spin 1/2, while the term, 'holon' refers to the bo- 
son of positive charge +e with spin 0. We find that the 
s-wave holon pairing ((&J&I)) but not the (i-wave holon 
pairing is stable for the d-wave spinon pairing ({fitfj\}), 
thus allowing the condensation of d-wave hole pairs in 
the language of 'hole'. In addition, it is shown that the 
spin gap remains nearly unchanged below the boson pair 
condensation temperature. 
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